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ABSTRACT 

The most direct signatures of particle acceleration in flares are 
energetic particles detected in interplanetary space and in the Earth’s 
atmosphere, and gamma rays, neutrons, hard X-rays, and radio 
emissions produced by the energetic particles in the solar at- 
mosphere. We review the stochastic and shock acceleration theories 
in flares, and we discuss the implications of observations on parti- 
cle energy spectra, particle confinement and escape, multiple ac- 
celeration phases, particle anistropies, and solar atmospheric 
abundances. 

1. INTRODUCTION 

Acceleration of energetic particles is a widespread phenomenon in nature, one 
that occurs at a variety of sites ranging from the Earth’s magnetosphere to 
distant objects such as supernovae, active galaxies, and quasars. There are, 
in fact, many explosive phenomena in astrophysics, solar flares among them, 
in which energetic particles are routinely produced and often contain a large 
fraction of all the available energy. 


47 

PRECEDING PAGE BLANK-NOT FILMED 



It is widely believed [e.g., Syrovatskii, 1981] that solar flares draw their energy 
from the annihilation of magnetic fields. The rapid energy deposition following 
the annihilation should be an important source of turbulence and shocks. As 
proposed by Fermi [1949], charged particles can be accelerated to high energies 
in repeated reflections from magnetized clouds, or, in the more recent view, 
from hydromagnetic turbulence and shocks. This mechanism must be impor- 
tant in solar flares where shocks are known to exist and turbulence is expected 
to be produced by both shocks and other mechanisms. In addition, particles 
may also be accelerated in the strong electric fields which should accompany 
the annihilation of magnetic fields. 

Indeed, there is ample evidence for particle acceleration in flares. Energetic 
particles accelerated at the Sun are directly detected in interplanetary space 
and the neutral radiations (radio emissions, hard X-rays, gamma rays, and 
neutrons) produced by the accelerated particles in the solar atmosphere are 
observed with detectors on the ground and in space. These observations pro- 
vide a variety of information on the properties of the particles, including their 
energy spectrum, total numbers at the Sun and in interplanetary space, and 
angular distribution at the Sun. In addition, the measured composition of 
solar flare particles in interplanetary space [e.g., Breneman and Stone, 1985] 
and flare gamma ray line spectroscopy [Murphy et al., 1985] provide new in- 
formation on the chemical composition of the solar atmosphere. 

In a previous paper [Forman, Ramaty, and Zweibel, 1986, hereafter FRZ] 
we have reviewed much of the data, as well as the stochastic, shock, and elec- 
tric field theories of particle acceleration in solar flares. Several other recent 
reviews are also available, e.g., Goldman and Smith [1986] on radio emis- 
sions, Dennis [1985] on hard X-ray bursts, Hudson [1985] on ions in solar 
flares, and Lin [1985] on solar electrons in the interplanetary medium. Here 
we discuss two of the leading acceleration mechanisms, the stochastic Fermi 
process and diffusive shock acceleration. These mechanisms not only are ex- 
pected to be important for particle acceleration in flares but also appear to 
be capable of accounting for many observed solar flare phenomena related 
to the acceleration of protons and relativistic electrons [Ellison and Ramaty, 
1985; Murphy, Dermer and Ramaty, 1987]. We also discuss recent observa- 
tions of charged particles, gamma rays, and neutrons, with particular focus 
on the simultaneous observations of flares in all three of these channels. 
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The present paper is the expanded version of the talk given by one of us (MAF) 
at Frank McDonald’s sixtieth birthday symposium. Frank, through his ex- 
tensive and pioneering observations of accelerated particles of planetary, solar, 
interplanetary, and galactic origins, has made fundamental contributions to 
the understanding of astrophysical particle acceleration. Frank’s most recent 
contribution to solar flare studies, the observation of accelerated particles from 
two neutron producing flares [McDonald and Van Hollebeke, 1985], is fun- 
damental to the considerations developed and discussed in the present paper. 


2. ACCELERATION MECHANISMS 
A. Stochastic Acceleration 

Processes in turbulent plasmas which cause particles to change their energy 
in a random way with many increases and decreases in energy lead to stochastic 
acceleration. In the original stochastic Fermi mechanism [Fermi, 1949], the 
process was idealized as reflection from randomly moving magnetized clouds. 
Stochastic acceleration can also result from resonant pitch-angle scattering 
from Alfven waves with wavelengths of the order of the particle gyroradius. 
To accelerate particles these waves must propagate both parallel and an- 
tiparallel to the average magnetic field [Skilling, 1975]. Other modes of 
stochastic acceleration, called magnetic pumping and transit-time damping, 
occur through interaction with magnetosonic waves whose wavelengths are 
much longer than the particle gyroradius [Kulsrud and Ferrari, 1971; Melrose, 
1980; and Achterberg, 1981]. These modes require additional pitch-angle scat- 
tering to keep the particles isotropic. Langmuir (plasma) waves or other elec- 
trostatic waves with phase velocities of the order of the particle speed will 
also accelerate particles stochastically [Melrose, 1980]. 

When the random energy increments are small compared to the particle energy, 
stochastic acceleration can be described as diffusion in momentum space 
characterized by a momentum diffusion coefficient D . Expressions for D_ 
for the various processes mentioned above were summarized in FRZ. 
Stochastic-acceleration spectra can be obtained from these D ’s by solving 
a transport equation [see FRZ and Droge and Schlickeiser, 1986] which takes 


49 


into account the injection of the particles, diffusion in momentum space, non- 
diffusive energy changes (e.g., ionization losses, and nuclear collisions), and 
escape from the acceleration region. Several of the published stochastic- 
acceleration spectra are given in FRZ. 

Here we discuss the simple model of acceleration by hard sphere scattering 
which has been applied extensively to both solar flare particle observations 
in interplanetary space [McGuire and von Rosenvinge, 1984] and gamma ray 
and neutron production models at the Sun [Murphy and Ramaty, 1984; Mur- 
phy, Dermer, and Ramaty, 1987]. In this treatment [Ramaty, 1979; FRZ] the 
scattering mean free path is assumed to be independent of particle energy and 
species, the acceleration region is characterized by a constant escape time and 
all additional losses are ignored. With a steady source of q particles cm -3 
s~‘ at an injection momentum p , the steady state particle density in phase 
space as a function of particle momentum p is given by: 


f = 6 q / (4 tt p 0 2 me a) (p Q /p) • 

• I 2 (2(3p 0 /(mc aT)) 1/2 ) K 2 (2(3p/(mc aT)) 1/2 ), (1) 


if both p and p 0 are nonrelativistic, and by: 

f = q / (4 tt p 0 3 a (1 +4/(3 aT)) ,/2 ) . 

,p/ Po )-3 /2 -V(974Twr) _ (2) 

when both p and p 0 are ultrarelativistic. Here m is particle mass, p particle 
momentum, T the escape time from the acceleration region, and a = (5V) 2 / 
(X c) where (5V) 2 is the mean square velocity of the scatterers and X the dif- 
fusion mean free path in momentum space. The combination of parameters 
aT characterizes the shape of the spectrum such that a larger value of aT 
corresponds to a harder spectrum. The phase space density f is related to the 
differential particle density, N(E), or the differential particle flux, J(E), N(E) 
= J(E)/v = A p 2 f/v, where E is energy per nucleon, A is particle mass 
number, and v particle velocity. In both Equations (1) and (2) p must be greater 
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than p 0 . For p < p 0 , the arguments of I 2 and K 2 in Equation (1) are inter- 
changed, and the minus sign in front of the square root in the exponent in 
Equation (2) is changed to a plus sign. 

The nonrelativistic expression, Equation (1), can be used to fit solar flare ac- 
celerated ion spectra up to about 100 MeV/nucleon. Such fits are discussed 
in Section 3. The ultrarelativistic expression. Equation (2), could be used to 
fit electron spectra above about 1 MeV. But if aT is rigidity independent, 
the values of aT that fit the proton spectra produce ultrarelativistic electron 
spectra which are much steeper than those observed [e.g., Ramaty, 1979]. 
Clearly, a complete stochastic acceleration model should take into account 
the rigidity dependence of aT, but no attempts have yet been made (see FRZ) 
to compare the theory with data in detail. 

A very important question in all particle acceleration theories, including 
stochastic acceleration, is that of injection. We first note that the basic con- 
cept of stochastic acceleration assumes that the energy changes are small com- 
pared with the particle energy and therefore the particle velocity must be much 
greater than 5V. Furthermore, for resonant scattering, ions must have v > 
V A to scatter from Alfven waves and electrons must have v > 43V A to scat- 
ter from whistlers [Melrose, 1974]. An additional injection condition is set 
by the requirement that the systematic acceleration rate due to diffusion in 
momentum space be larger than the ionization and Coulomb energy loss rates 
of the particles. A detailed comparison of these rates was given in FRZ. 

The acceleration time in stochastic acceleration can be studied from the time 
dependent solutions of the transport equation. The relevant parameters are 
5V and X . If we set SV = V A , then a typical value is 5V = 1000 km/sec, 
corresponding to a magnetic field of 100 gauss and an electron density of 
5 x 10 10 cm -3 . The diffusion mean free path should be at least as large as the 
particles gyroradius. Taking X to be 100 times the gyroradius of a 20 MeV 
proton, we obtain X = 7 x 10 5 cm. Then the acceleration parameter a = 0.5. 
For this a, the time dependent solutions of the transport equation obtained 
[Ramaty, 1979; FRZ] for impulsive injection and perfect trapping (T— oo) 
imply that a significant number of protons are accelerated to tens of MeV 
in less than 1 sec. In general, a time of the order 1/a is required to accelerate 
particles to relativistic energies. 
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B. Shock Acceleration 


Solar flare shocks propagate upward through the solar corona at speeds of 
about 500 to 2000 km/sec, as indicated by Type II radio bursts [e.g., Goldman 
and Smith, 1986] and laterally through the chromosphere where they are seen 
as Moreton waves [Uchida, 1974]. The occurrence of solar energetic particles 
in space is strongly correlated with flares having Type II bursts [Svestka and 
Fritzova, 1974], Flare shocks, corotating shocks, and planetary bow shocks 
are observed to accelerate particles in interplanetary space (see references in 
FRZ). A flare shock can transport particles in an energy independent man- 
ner through the corona until they escape onto open field lines. Shock accelera- 
tion has been reviewed by Toptygin [1980] and Axford [1982] and applied 
to solar flares by Achterberg and Norman [1980], Decker, Pesses, and Arm- 
strong [1981], Ellison and Ramaty [1985], and Lee and Ryan [1986]. 

There are basically two types of shock acceleration: scatter-free, in which par- 
ticles gain energy by reflection in a single shock encounter [Sonnerup, 1973; 
Armstrong, Pesses, and Decker, 1985; references in FRZ] and diffusive, in 
which particles gain energy by repeated scattering between the converging 
plasmas on either side of the shock [e.g., Axford, Leer, and Skadron, 1977; 
Axford, 1982 and; Forman and Webb, 1985]. The scatter-free mechanism can 
enhance the particle energy by about an order of magnitude if the shock is 
nearly perpendicular (i.e. , the magnetic field is nearly perpendicular to the 
shock normal), but in that case only particles with speeds which are already 
much greater than the shock speed can be reflected. Further acceleration, 
however, requires multiple reflections. These are possible if there is particle 
scattering in the fluid flow or if the particles are trapped between converging 
shocks in a flare loop [Wentzel, 1965; Bai et al., 1983]. 

Acceleration by diffusive scattering across the shock is a first order Fermi 
process, in the sense that every shock crossing results in an energy gain. It 
is in principle more efficient than stochastic acceleration because it derives 
energy directly from the compression of the flow at the shock. For this 
mechanism to be effective, there must be adequate particle scattering both 
upstream and downstream of the shock. The passage of the shock is expected 
to generate turbulence downstream which will scatter the particles. Scatter- 
ing upstream, however, is more problematic [Holman, Ionson, and Scott, 
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1979]. Observations [Tsurutani and Rodriquez, 1981] of interplanetary shocks 
and planetary bow shocks show that when they are nearly parallel there is 
a very turbulent foreshock region capable of scattering particles. Such a region 
could be produced by the accelerated particles themselves [Achterberg and 
Norman, 1980; Lee, 1982]. 


In the simplest example of diffusive shock acceleration at a planar shock where 
the only losses are due to convection of the particles away from the shock 
downstream, the energetic particle density in phase space is given by a power 
law in momentum, f ~ p~", where a = 3V/AV, V is the shock speed and 
AV the discontinuity in the plasma speed at the shock. In terms of the shock 
compression ratio, r, a = 3r/'(r-l). For a strong shock in a nonrelativistic 
fluid r = 4 and hence f ~ p~ 4 . For weaker shocks, 4 > r > 1 and the power 
law is steeper. The momentum power law p~" implies [Ellison and Ramaty, 
1985] that differential flux as a function of energy per nucleon, J(E), has an 
energy dependent spectral index which approaches s = (r + 2)/(r- 1) in the 
ultrarelativistic regime and s/2 in the nonrelativistic regime. 


A variety of effects truncate the power law behavior of shock-accelerated spec- 
tra at high energies. The effects that have been considered specifically are 
adiabatic deceleration [Lee and Ryan, 1986], shock lifetimes comparable to 
particle acceleration times [Forman, 1981] and shock sizes comparable to par- 
ticle diffusion lengths [e.g., Ellison, 1984]. These effects all produce spectra 
that turn over at high energies when the diffusion coefficient increases with 
momentum. The derivation of the exact forms of these turnovers generally 
requires numerical treatments. In many cases, however, they can be approx- 
imated by exponentials [e.g., Ellison and Ramaty, 1985], 


E(P) ~ P~ a exp(-E/E 0 .), 


( 3 ) 


where f. is the phase space density of particles of species i, and E^ is the turn- 
over kinetic energy for species i. In the case of particle escape at an escape 
boundary at a given distance from the shock, E M depends on the distance to 
this boundary and the diffusion coefficent. If the diffusion mean free path 
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is proportional to particle gyroradius, the turnover energies of the various 
species are related by: 


v(E 0i ) R(E„) = v(E 0p ) R(E 0p ), (4) 

where R is particle rigidity, E is the proton turnover kinetic energy, and 
is energy per nucleon for nuclei and energy for electrons. Equation (3) can 
be used to fit both ion and electron spectra. We present such a fit in the next 
section. 

As for stochastic acceleration, and for shock acceleration as well, the ques- 
tion of injection is very important. Ionization and Coulomb energy losses to 
the ambient medium have the same role in determining injection conditions 
in shock acceleration as they do in stochastic acceleration. In addition, for 
diffusive shock acceleration, particles downstream must have sufficient velocity 
to overtake the shock. This is at least (V-AV) directed towards the shock, 
and increases as (cos ^ )' 2 , where 4J is the angle between the downstream 
field and the shock normal. The velocity V - AV is at least as great as V A or 
AV, and with the additional (cos 4 1 )~ 2 factor, the threshold for shock ac- 
celeration could be higher than for stochastic acceleration. 

Another injection condition is set by the finite width of the shock which could 
depend on many parameters including the pressure of the accelerated par- 
ticles. When this pressure is taken into account [Axford, Leer, and Skadron, 
1977; references in FRZ] all or part of the velocity change AV is smoothed 
out over a length scale ~ k/V, where k is the diffusion coefficient of the 
energetic particles averaged over their energy spectrum. This smoothing could 
affect the composition of the accelerated particles [Eichler, 1979]. Drury, Ax- 
ford, and Summers [1982] show analytically that when k is independent of 
energy, the smoothing causes the dominant accelerated species (i.e., protons) 
to have a steeper spectrum than in the case of an infinitely thin shock. Minor 
species which are only partially stripped could have larger diffusion coeffi- 
cients than the protons if the diffusion mean free path is rigidity dependent, 
and therefore for them k > k. Drury, Axford, and Summers [1982] show that 
the spectrum of such minor species is flatter than for protons and approaches 
that of an infinitely thin shock. 
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The mean rate of energy gain in shock acceleration is (see FRZ): 

dE/dt = 3/4 (V AV) (p/X). (5) 

Ellison and Ramaty [1985] have shown that if X is proportional to gyroradius, 
the acceleration time for both protons and electrons is: 

t a = 2x10 “ f (E-E. nj )/((B/100 gauss)(V AV/10 6 km 2 /sec 2 )), (6) 

where t a is in sec, E and E. nj are the final and injection energies in MeV, and 
f is the ratio of X to the gyroradius. It is expected that f is much greater than 
unity. For f = 100, B = 100 gauss and V = AY = 10 3 km/sec, the accelera- 
tion time to 20 MeV is ~ 0.4 sec. We note that the acceleration time is the 
same for protons and electrons, a result which follows from the assumed linear 
dependence of the mean free path on gyroradius. 


3. OBSERVATIONS AND THEIR INTERPRETATIONS 

We discuss the following: (1) energy spectra, (2) the trapping and escape of 
the particles and the related problem of multiple acceleration phases, (3) 
angular distributions at the Sun, and (4) abundance determinations from 
energetic particle and gamma ray observations. 

A. Energy Spectra 

The energy spectrum of accelerated particles is perhaps the most important 
diagnostic of acceleration mechanisms. Information on energy spectra is ob- 
tained from both the charged particles observed in interplanetary space, and 
the gamma ray and neutron data. The interplanetary observations provide 
information on the energy spectrum of the particles which escape from the Sun. 

As indicated by upper limits on 2 H, 3 H, Li, Be, and B [McGuire, von Rosen- 
vinge, and McDonald, 1979], the amount of matter traversed by the escaping 
particles (< 0.1 g/cm 3 ) is generally insufficient for the production of the 
observed gamma rays [Ramaty, 1986]. Gamma rays and neutrons are pro- 
duced by the particles which slow down and thermalize in the solar atmosphere. 
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Observations of these neutral emissions, therefore, provide information on 
the particles which remain trapped at the Sun. 

The energy spectra obtained from charged particle observations have been 
reviewed by McGuire and von Rosen vinge [1984]. In determining spectra from 
such observations the effects of coronal and interplanetary transport must 
be minimized. It was shown [Van Hollebeke, MaSung, and McDonald, 1975] 
that this can be achieved by considering only particle events from flares that 
are well connected magnetically to the observing spacecraft and by construc- 
ting the energy spectra at times of maximum intensity at each energy. Using 
this technique, McGuire, von Rosenvinge, and McDonald [1981] and McGuire 
and von Rosenvinge [1984] showed that over the broad energy range (from 
about 1 to 400 MeV) proton spectra can be fit with the Bessel function spec- 
tra expected from stochastic acceleration [Equation (1)] with values of aT 
between 0.015 and 0.035. But as emphasized in FRZ, energy spectra obtained 
from shock acceleration could also fit these data. 

Proton spectra have been determined from gamma ray and neutron observa- 
tions using three different techniques [Murphy and Ramaty, 1984; Murphy, 
Dermer, and Ramaty, 1987]. The first one considers the ratio of the fluence 
in excess of a power law in the 4-7 MeV range to the fluence in the 2.223 
MeV line. This excess is due to nuclear line emissions. The second technique 
uses observed neutron spectra which depend on the spectrum of the protons 
which produce the neutrons; and the third is based on the ratio of the 4-7 
nuclear excess to the pion decay emission at photon energies > 10 MeV. The 
first technique is sensitive to the spectrum in the 10-100 MeV range, the sec- 
ond probes the spectrum in the 50 MeV - several GeV range and the third 
is sensitive to the ratio of the number of particles in the 10 - 100 MeV range 
to the number above the pion production threshold (few hundreds MeV). 

When analyzed in terms of the Bessel function of Equation (1), the 4-7 MeV 
to the 2.223 MeV line ratios for 15 flares yielded values of aT in the range 
from 0.014 to 0.038 [Murphy and Ramaty, 1984; Hua and Lingenfelter, 1987a]. 
This range is essentially the same as that obtained from the interplanetary 
charged particle data. However, the implications of this overlap are not clear 
because most of the particle events correspond to a different set of flares than 
the one used to derive the gamma ray results. 
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There are, in fact, only three flares for which energy spectra have been deter- 
mined simultaneously from gamma ray and charged particle data: June 7, 
1980; June 21, 1980; and June 3, 1982. For the June 7 flare, the 4-7 MeV 
to 2.223 MeV ratio yields aT = 0.021 ±0.003 [Murphy and Ramaty, 1984; 
Hua and Lingenfelter, 1987a] and the charged particle data yield aT = 0.013 
[McGuire and von Rosenvinge, 1985]. For the June 21 flare, the 4-7 MeV 
to 2.223 MeV ratio yields aT = 0.022 ± 0.007 [Hua and Lingenfelter, 1987a] 
and the charged particle data yield aT = 0.025 [Murphy, Dermer, and Ramaty, 
1987] using the data of McDonald and Van Hollebeke [1985] shown in Figure 
1. For the June 21 flare aT was also determined [Murphy and Ramaty, 1984] 
from neutron observations [Chupp et al., 1982] yielding aT = 0.025. We see 
that for both the June 7 and June 21 flares the spectrum of the particles which 
escaped from the Sun could have been similar to the spectrum of the par- 
ticles which slowed down and thermalized in solar atmosphere. It is possible 
that for these flares a common acceleration mechanism operating at a com- 
mon acceleration site is responsible for both particle populations. 


The interplanetary proton and electron spectra from the June 3, 1982 flare 
[McDonald and Van Hollebeke, 1985] are shown in Figure 2 [from Murphy, 
Dermer, and Ramaty, 1987]. As can be seen, a power law fits the proton spec- 
trum better than does a Bessel function. This power law is given by Equation 
(3) with s = 2.4 (a =4.4) and E 0p »300 MeV. The gamma ray data from 
this flare, however, cannot be reproduced with a single, time independent, 
proton spectrum [Murphy, Dermer, and Ramaty, 1987]. Early in the flare 
(for about 100 sec) the proton spectrum can be fit with a Bessel function with 
aT = 0.04, but later in the event the observed ratio of the 4-7 MeV nuclear 
excess [Prince et al., 1983; Chupp et al., 1987] to the emission from the decay 
of pions [Forrest et al., 1985, 1987] implies that the spectrum must harden. 
At these later times the spectrum is consistent with the proton spectrum ob- 
served in interplanetary space (s = 2.4 and E 0p > 300 MeV). 

The inferred hardening of the proton spectrum suggests that not all the par- 
ticles accelerated in the June 3 flare had a common origin. Murphy, Dermer, 
and Ramaty [1987] show that a wide variety of data from the June 3 flare 
(pion decay emission, nuclear deexcitation lines, the 2.223 MeV and 0.511 
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Figure 1. Interplanetary proton and electron observations [McDonald and 
Van Hollebeke, 1985] of the June 21, 1980 flare. The solid curve corresponds 
to a stochastic-acceleration spectrum with aT = 0.025, and the dashed curves 
correspond to shock-acceleration spectra with s = 3.3 and E 0p -30 MeV for 
the protons and E oe = 59 MeV for the electrons. 
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Figure 2. Interplanetary proton and electron observations [McDonald and 
Van Hollebeke, 1985] of the June 3, 1982 flare. The solid curve corresponds 
to a stochastic-acceleration spectrum with aT= 0.04, and the dashed curves 
correspond to shock-acceleration spectra with s=2.4 and E 0p and E w tending 
to infinity. 
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MeV lines, neutrons and interplanetary charged particles) can be explained 
by a model which incorporates two distinct particle populations with different 
acceleration histories and different but time independent energy spectra. They 
suggest that the softer proton spectrum is accelerated early in the flare by 
the first of two acceleration phases, while later in the flare second phase acceler- 
ation produces a much harder spectrum. We shall discuss these multiple ac- 
celeration phases in the next subsection, but before that we will comment on 
the maximum energy of particles accelerated in flares. 

As just mentioned, the interplanetary proton spectrum of the June 3, 1982 
flare was exceptionally hard, but these observations determined the proton 
spectrum only up to ~200 MeV. The ground-based neutron observations from 
this flare [Debrunner et al., 1983], which are very sensitive to protons of GeV 
energies, set an upper limit, E 0p < 3 GeV [Murphy, Dermer, and Ramaty, 
1987]. The highest energy particles from flares have been seen with ground- 
based detectors. These show that the interplanetary proton spectrum from 
solar flares occasionally extends to energies > 10 GeV [Debrunner et al., 1984], 


B. Trapping and Escape of Particles and Multiple Acceleration Phases 

With a few exceptions (the flares of August 4, 1972 and June 3, 1982), the 
number of protons deduced from gamma ray and neutron observations ex- 
ceeds by at least an order of magnitude the number obtained from in- 
terplanetary observations of the same flares. This implies that the gamma rays 
and neutrons are produced by particles accelerated in and confined to closed 
magnetic configurations, probably magnetic loops. On the other hand, there 
are many flares which produce large fluxes of interplanetary particles without 
producing detectable gamma rays or neutrons. The flare of December 9, 1981 
is an example [Oliver et al., 1983]. The number of escaping particles from 
this flare exceeds the upper limit set [Vlahos et al., 1987] by gamma ray obser- 
vation by a factor of 5. These particles are most likely accelerated at coronal 
sites with ready access to interplanetary space. 

These results are consistent with the two phase acceleration model suggested 
originally by Wild, Smerd, and Weiss [1963]. Recently Cane, McGuire, and 
von Rosen vinge [1986] have related the two acceleration phases to two classes 
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of solar flare particle events which can be distinguished by a variety of 
signatures including the duration of their associated soft X-ray emission 
[Pallavicini, Serio, and Vaiana, 1977]. The first class, characterized by short 
duration impulsive events, is probably due to impulsive energy release in com- 
pact source regions relatively low in the corona. The second class, characterized 
by longer duration soft X-ray emission, corresponds to energy release in ex- 
tended regions high in the corona. It is reasonable to associate the bulk of 
the gamma ray flares (including the June 7, 1980 and June 21, 1980 flares 
for which only a small fraction of the protons escaped from the Sun) with 
the first class or the first phase, and the interplanetary proton events which 
do not show detectable gamma rays (e.g., the December 9, 1981 flare) with 
the second class or phase. A similar classification has also been proposed re- 
cently by Bai [1986]. 


Of all flares, only the June 3, 1982 flare produced observable gamma ray emis- 
sion which can be associated with second phase acceleration. As in other gam- 
ma ray flares, this flare also produced nuclear line emission and neutron cap- 
ture radiation [Prince et al., 1983; Chupp et al., 1987], These emissions are 
typical of first phase acceleration [Murphy, Dermer, and Ramaty, 1987]. But 
the June 3 flare also produced emission from the decay of pions [Forrest et 
al., 1987] and the time profile of the observed pion decay emission strongly 
suggests second phase acceleration [Ramaty, Murphy, and Dermer, 1987]. 
The number of escaping particles from this flare was smaller by about an order 
of magnitude than the number deduced from the gamma ray and neutron 
data in the first phase and exceeded by a factor of about 50 the number deduced 
from the pion decay gamma rays in the second phase. In addition, as dis- 
cussed above, the spectrum of the escaping particles is much more consistent 
with the spectrum of the second phase particles than with that of the first 
phase ones. These results suggest that during its first phase the June 3 flare 
was similar to the other gamma ray flares in which the bulk of the interacting 
particles remained trapped and thermalized at the Sun. During the second 
phase of this flare, the bulk of the particles escaped. This phase of the June 
3 flare resembles the flares which produce interplanetary particles without 
detectable gamma rays and neutrons. Second phase gamma rays were seen 
from the June 3 flare only because of its very hard proton spectrum which 
led to efficient pion production. 
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C. Energetic Particle Anisotropies in the Solar Atmosphere 

The observations that continuum gamma ray emission is seen preferentially 
from flares close to the solar limb [Rieger et al., 1983; Vestrand et al., 1987] 
suggest that the angular distribution of the relativistic electrons in the region 
where the gamma rays are produced is anisotropic. It was suggested that this 
limb brightening could be due to an electron distribution that at injection was 
isotropic in the downward hemisphere [towards the photosphere, Petrosian, 
1985], and by electron distributions in the interaction region which were peaked 
at directions either parallel or perpendicular (downward) to the photosphere 
[Dermer and Ramaty, 1986]. 


The angular distribution of the protons in the gamma ray production region 
can be studied by comparing the 2.223 MeV line, which is produced by neutrons 
moving downward into the photosphere, with the neutron flux observed near 
Earth. A recent study [Hua and Lingenfelter, 1987b] indicates that both an 
isotropic ( 4 t ) proton distribution and a distribution peaking at directions 
parallel to the photosphere are consistent with the data, but a distribution 
which is isotropic in only the downward hemisphere is not. This result, coupled 
with those on relativistic electrons, suggests that the gamma rays are probably 
produced by a mirroring distribution peaking at directions parallel to the 
photosphere. 


D. Solar Elemental Abundances 


The observation of energetic particles escaping from the Sun and gamma ray 
lines from solar flares have provided two new techniques for determining 
elemental abundances in the solar atmosphere. The elemental composition 
of the energetic particle depends on the composition of the ambient medium 
from which these particles are accelerated, but the injection and acceleration 
processes are expected to modify the composition significantly. Gamma ray 
line intensities are directly proportional to the abundance of elements in the 
ambient solar atmosphere, but so far the gamma ray spectrum of only one 
flare was analyzed in sufficient detail to determine abundances. 
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Nuclei heavier than He in solar energetic particles were first detected by Fichtel 
and Guss [1961] and since then many measurements of such particles have 
been made (see references in FRZ). The observed solar energetic particle com- 
position is highly variable. The most dramatic departure of a solar energetic 
particle abundance from its photospheric value is that of 3 He [Garrard, 
Stone, and Vogt, 1973]. Here very large enhancements are observed in the 
3 He/ 4 He ratio above its likely photospheric value [see Kocharov and 
Kocharov, 1984 for review]. These enhancements are probably caused by reso- 
nant heating of 3 He in the ambient atmosphere which could lead to the 
preferential accelerations of 3 He [Fisk, 1978; Kocharov and Kocharov, 
1978]. It has been shown recently [Reames, von Rosenvinge, and Lin, 1985] 
that the preferential acceleration of 3 He is most likely a first phase 
phenomenon. 

Observations of the 2.223 MeV line from solar flares have been used to deter- 
mine the photospheric abundance of 3 He. It was shown by Wang and Rama- 
ty [1974] that the photospheric 3 He is an important neutron sink and 
therefore the observed 2.223 MeV line, which results from neutron capture 
on 'H, can set limits on the photospheric 3 He abundance. Recently, Hua and 
Lingenfelter [1987c] used the detailed June 3, 1982 observations of the 2.223 
MeV line [Prince et al., 1983] to determine the 3 He abundance. They 
obtained 3 He/H = (2.3 ±1.2) x 10" 5 . 

It would appear that the large variations of the observed composition of the 
solar energetic solar particles would preclude the determination of abundances 
in the ambient medium. Nonetheless, it has been suggested [Meyer, 1985a] 
that for non 3 He rich flares (for which the abundance variations are rela- 
tively small) the acceleration and injection effects could be separated from 
effects related to the ambient medium composition. In particular, it has been 
shown [Breneman and Stone, 1985] that the ratio of the abundance of elements 
observed in individual flares to the mean abundance (determined by averag- 
ing abundances of several flares) is a monotonic function, (Q/M) x , where Q 
and M are ionic charge and mass, and x varies from flare to flare. These mean 
abundances could be similar to those of the ambient medium. 

In Table 1 we show abundances derived from various sources. The local galactic 
(LG) set [Meyer, 1985b] is thought to represent photospheric abundances, 
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TABLE 1 


Elemental Abundances 

1 2 3 


Element 

Local 

1 

Corona 


SEP 


Gamma Rays 


Galactic 

1 








C 

1260.00 

(1, 

.26) 

600.00 

(3 

00) 

271 .00+- 

26.00 

288.00+ 

50.00 

N 

225.00 

(1. 

.41) 

100.00 

(1 

70) 

77.50+- 

5.20 

117.00+ 

91.00 

0 

2250.00 

(1. 

.25) 

630.00 

(1 

60) 

623. 00+ - 

35.00 

422.00+ 

62.00 

Ne 

325.00 

(1 

.50) 

90.00 

(1 

60) 

88 . 70+ - 

8.70 

199.00+ 

27.00 

Na 

5.50 

(1. 

.18) 

7.00 

(1 

70) 

7.33+- 

0.70 

- - 

- 

Mg 

105.00 

(1, 

.03) 

95.00 

(1 

30) 

120.60+- 

6.20 

68.00+ 

25.00 

Al 

8.40 

(1, 

.05) 

7.00 

(1 

70) 

8.74+- 

0.42 

- - 

- 

Si 

100.00 

(1 

.03) 

100.00 

<1 

30) 

100.00 


100.00+ 

28.00 

P 

0.94 

(1 

.24) 

- - 

- 


0.46+- 

0.07 

- - 

- 

S 

43.00 

(1 

.36) 

22.00 

(1 

70) 

22.20+- 

0.75 

48.00+ 

83.00 

Cl 

0.47 

(1 

.60) 

• • 

- 


0.21+- 

0.07 

- - 

- 

Ar 

10.70 

(1 

.50) 

5.40 

(1 

80) 

2.07+- 

0.33 

- - 

- 

K 

0.34 

(1 

.14) 

- - 

- 


0.33+- 

0.15 

- - 

- 

Ca 

6.20 

(1 

.14) 

7.50 

(1 

60) 

6.80+- 

1.10 

17.00+ 

15.00 

Ti- 

0.27 

(1 

.16) 

• - 

- 


0.38+- 

0.11 

- - 

- 

er 

1.29 

(1 

.10) 

• - 

- 


1 .43+- 

0.27 



Mn 

0.77 

(1 

.24) 

• - 

- 


0.52+- 

0.25 



Fe 

88.00 

(1 

■ 07) 

100.00 

(1 

50) 

95.90+- 

10.00 

76.00+- 

■18.00 

Ni 

4.80 

(1, 

.13) 

5.50 

(1 

70) 

3.38+- 

0.50 

. . . 

■ • 

Zn 

0.10 

(1. 

.22) 

- - 

- ■ 


0.11+- 

0.05 

. . . 

. . 


1. Meyer (1985b) 

2. Breneman and Stone (1985) 

3. Murphy et al. (1985) 
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even though the Ne abundance has not yet been measured in the photosphere. 
The coronal (COR) set, based on spectroscopic observations of the corona, 
is also from Meyer [1985b]. The mean solar energetic particle abundances 
(SEP) are from Breneman and Stone [1985]. The gamma ray (GR) set is from 
Murphy et al. [1985]. 

To compare these abundances, we have renormalized each pair of sets using 
multiplicative factors determined by minimizing x 2 - The ratios of elemental 
abundances for each pair are shown in Table 2; together with the number 
of degrees of freedom, n; the value of the reduced x 2 . x„ 2 ; and the corre- 
sponding probability. Here the elements are ordered by decreasing first ioniza- 
tion potential. As can be seen, it is very improbable that the LG and SEP 
sets are drawn from the same underlying population. In particular, the SEP 
to LG ratios for elements with high ionization potential are lower than the 
ratios for elements with low ionization potential. A similar suppression has 
been found in the coronal abundances relative to the local galactic abundances 
(see the COR/LG ratios), although here the difference is much less signifi- 
cant. As we have argued above, and as suggested by charge state measurements 
[Gloeckler, 1985], the solar energetic particles are probably accelerated in the 
corona. The SEP abundance set could therefore represent coronal abundances. 
The differences between coronal and photospheric abundances could be caused 
by charge dependent mass transport from the photosphere to the corona. Since 
the photosphere is collisionally ionized at a relatively low temperature, the 
transport could be less efficient for elements with high ionization potential, 
leading to suppressed abundances for such elements. 

The gamma ray set probably represents chromospheric abundances in a flare 
loop [Murphy and Ramaty, 1984]. As for the SEP/LG ratios, the GR/LG 
ratios for O and C are lower than the ratios for Mg, Si, and Fe. The ioniza- 
tion potentials of C and O are higher than those of Mg, Si, and Fe. The sup- 
pression of the C and O abundances could also be caused by charge depen- 
dent transport, from the photosphere to the chromosphere in this case. 
However, if the Ne abundance in the photosphere (where it cannot be 
measured) is the same as in the local galactic set, then the difference between 
the chromospheric and photospheric abundances must be due to additional 
processes, because correlation with first ionization potential alone would 
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TABLE 2 


Ratios of elemental abundances for the six combinations of the data sets of Table 1. 

2 

Also shown are the degrees of freedom, n, the minimal reduced X and the corresponding 
probability, P(X), for each set of ratios. The elements are ordered by decreasing first 
ionization potentials. 


Element SEP/LG 

GR/LG 

SEP/COR 

GR/COR 

GR/SEP 


COR/ LG 

Ne 

0.27+-0.14 

0.80+-0.41 

0.93+ -0.57 

1 .84+- 1 . 13 

2.17+-0. 

36 

0.32+-0.25 

Ar 

0. 19+-0.10 

.... 

0.36+ -0.30 

.... 



0.58+ -0.55 

N 

0.34+-0.14 

0.68+-0.59 

0.73+-0.52 

0.97+- 1 .02 

1 .46+- 1 . 

14 

0.51+-0.41 

0 

0.27+0.07 

0.24+ -0.07 

0.93+-0.56 

0.56+-0.34 

0.66+-0. 

10 

0.32+-0.21 

Cl 

0.42+ -0.29 


.... 

.... 


- 

.... 

C 

0.21+-0.06 

0.30+-0.09 

0.43+ -0.85 

0.40+ -0.80 

1.03+-0. 

20 

0 . 55+ -1.10 

P 

0.48+-0. 13 

.... 

.... 

.... 



.... 

s 

0.50+-0. 18 

1 .45+2.56 

0.95+-0.67 

1.81+-3.39 

2.09+-3. 

62 

0.59+-0 .46 

Zn 

1.05+-0.57 




. . . . 



Si 

0 . 97+ -0.03 

1.30+-0.37 

0.95+-0.28 

0.83+-0.34 

0.97+-0. 

27 

1 . 15+-0.35 

Fe 

1.06+-0.13 

1.12+-0.28 

0.91+-0.46 

0.63+-0.35 

o 

£ 

o 

20 

1.30+-0.66 

Mg 

1.12+0.07 

0.84+-0.31 

1.20+-0.37 

0.60+-0.28 

0.55+-0. 

20 

1 .04+ -0.31 

Ni 

0.69+-0.13 

.... 

0.58+-0.42 

.... 



1 .32+0.94 

Mn 

0.66+-0.35 



.... 

.... 

- 

.... 

Ti 

1.37+-0.45 

* - - - 


.... 



.... 

Cr 

1 . 08+ *0.23 

.... 

.... 

.... 



.... 

Ca 

1.07+-0.23 

3.56+-3.18 

0. 86+ - 0.53 

1 .89+- 2.01 

2.42+*2. 

.17 

1 .39+- 0.86 

Al 

1.01+-0.07 

.... 

1.18+-0.83 

.... 


- 

0.96+ -0.67 

Na 

1.30+-0.26 

.... 

0.99+-0.70 

.... 



1 .46+- 1 .06 

K 

0.94+-0.45 

.... 

.... 


. . . 

. 

.... 


n 

19 

8 

12 

8 

8 

12 

2 

X (n) 

2.29 

2.47 

0.15 

0.51 

3.61 

0.85 

P 

-3 

1.1x10 

-2 

1.1x10 

>0.99 

0.85 

-3 

<10 

0.6 
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predict a lower Ne abundance than the one obtained from the gamma ray 
studies. 

Also shown in Table 2 are ratios of the gamma ray set to the SEP and co- 
ronal sets, and the ratios of the coronal set to the local galactic set. The dif- 
ference between the GR and SEP sets is significant. In addition to the dif- 
ferent Ne abundances, we also note that Mg in the GR set is significantly lower 
than in the SEP set (see also Table 1). There is no simple interpretation for 
this result. Clearly, the effects of injection and acceleration on the mean SEP 
abundances are not fully understood. Also, gamma ray data from more flares 
are needed to establish the constancy or variability of the abundances derived 
from gamma ray spectroscopy. 

Elemental compositions have also been studied for 3 He rich flares. Mason 
et al. [1986] found that the heavy ion enrichments seen in these flares are un- 
correlated with the 3 He enrichment. As pointed out above, the 3 He enrich- 
ment is most likely due to preferential heating and acceleration. Mason et 
al. [1986] suggest that the heavy ion enrichments seen in 3 He flares are 
caused by enrichments in the ambient coronal composition. 


4. SUMMARY 

That particle acceleration plays a dominant role in the physics of solar flares 
has been known for some time. In the present paper we have emphasized 
phenomena related to the acceleration of ions and relativistic electrons. We 
treated the presently best understood acceleration mechanisms, stochastic ac- 
celeration, and shock acceleration. We presented recent charged particle, gam- 
ma ray, and neutron observations, and we discussed their implications. These 
observations are highly complementary. The gamma rays and neutrons pro- 
vide information on the particles which thermalize in the solar atmosphere, 
while the charged particles observations provide a direct sample of the escap- 
ing particles. The combination of information on the escaping and trapped 
particle populations strongly suggests the existence of multiple acceleration 
phases, which could be the manifestation of multiple acceleration sites, times, 
and perhaps even mechanisms. 
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Particles accelerated in solar flares also provide information on elemental abun- 
dances in the solar atmosphere. This information is obtained by measuring 
the composition of the escaping particles, by observing the gamma ray lines 
which are produced from ambient nuclei excited by the accelerated particles, 
and by studying the 2.223 MeV line which gives information on the 3 He in 
the photosphere. Abundance variations in the chromosphere and corona are 
suggested by these observations. 
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